Strong effect of surfaces on resolution limit of negative-index "superlens 
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We show that subwavelength imaging by negative index materials (NIM), related to their "soft" 
electromagnetic response, is very (and non trivially) sensitive to the surface properties. A minute 
deviation of dielectric permittivity e or magnetic permeability /i from the ideal values e — fi = ~1 
in thin surface layer(s) results in drastic reduction of the resolution limit of a NIM slab. There may 
be a gap in the polariton spectrum and this would allow establishment of a stationary regime even 
without losses. 



Negative- index materials (NIMs), discussed theoreti- 
cally by Veselago in 1960s [1], with negative values of 
both the dielectric permittivity, e < 0, and permeability 
/i < 0, should show a negative refraction index in the 
Snell's law. This is a result of having the opposite phase 
and group velocities in NIM, as immediately follows from 
the Maxwell's equations. This behavior is very general 
and was noticed by Maldelshtam already in 1945 [2] (see 
discussion in [3]). NIMs have not been found in nature 
yet, but artificial metamaterials (metallic wire structures, 
practically, antenna arrays) have been recently shown 
to have a negative refraction in the microwave region 
[4]. The interest in metamaterials has intensified since 
Pendry's speculations about "superlensing" with the use 
of the NIMs [5]. 

The term "superlensing" , or subwavelength imaging, 
means a restoration of evanescent waves (exponentially 
decaying away from the source), that are getting ampli- 
fied inside NIM and, together with propagating waves, 
perfectly restore the image. However, this happens at 
extremely special conditions. To begin with, the perfect 
restoration of image takes place for ideal NIM only, i.e. 
for NIM with e = ^ = — 1. It is easy to see from geomet- 
rical optics [1] that a parallel slab of ideal NIM produces 
an exact replica of a source if the source is closer to the 
slab surface than the slab thickness. However, if e/i ^ 1 
this property does not exist [6], i.e. even within the geo- 
metrical optics the replica will not be perfect. We shall 
assume below that e/i = 1 , and will discriminate between 
the above defined ideal case and a non-ideal one where 
efj, = 1 but e, fj, —1. Although these two cases are 
equivalent within geometrical optics, they are essentially 
different within the wave optics. 

Full wave optics consideration of the radiation prop- 
agation through ideal NIM slab has been done in Ref. 
[5]. It has been shown that evanescent waves decaying 
while propagating from the source transform into increas- 
ing ("anti-evanescent") ones in the slab. Moreover, this 
amplification is such that the evanescent components are 
exactly restored in the replica. It could be said, therefore, 
that the replica remains perfect within the wave optics 



too. A difficulty with this result was that the stationary 
state supposed in [5] could not be achieved without ab- 
sorption [7] , and the problem should be reformulated as a 
non-stationary one [8]. Apart from this, a non-ideal slab 
with e/j, — 1 has been shown to not generate the perfect 
replica [9,10]. The case of homogeneous e,fi has been 
considered by Haldane [9]. We will argue in this letter 
that even much more delicate non-ideality, namely, when 
e, /i 7^ —1 in very thin surface layers only, and the thick- 
ness of the layers being much less that the wavelength, 
the replica will be distorted similarly to the case of the 
bulk non-ideality. We will also show that the surface 
non-ideality, unlike the bulk one, may change the char- 
acter of the amplification of the evanescent waves making 
possible a stationary state even without absorption. 

The physical reason of amplification of evanescent 
waves by NIM was nicely commented by Haldane [9] . The 
mechanism is valid not only for NIM's but for any slab 
with surfaces supporting surface electromagnetic ("po- 
lariton") waves, i.e. for any slab with negative e (/i) sup- 
porting polariton withp- (s-) polarization. The key point 
is that the incident radiation may be in resonance with a 
polariton mode of the NIM slab. Evanescent waves with 
wavevectors along the surface smaller than the wavevec- 
tor of the resonant polariton are amplified inside the slab, 
while the rest is not. The amplitude of the stationary os- 
cillations associated with the resonant mode will formally 
be infinite, and the contribution of the Fourier harmonics 
close to this mode will be exaggerated in the image. Con- 
sequently, the replica would be distorted. The ideal NIM 
is an exception, since there the resonance takes place at 
an infinite wave vector, and the exponentially large in- 
crease of the evanescent waves is fully compensated by an 
exponential decrease of the evanescent wave propagating 
from the slab surface to the replica image [9]. Moreover, 
the amplification of the (non-resonant) evanescent waves 
with finite wavevectors is such that they are exactly re- 
stored in the replica. Bearing in mind that, in addition, 
the propagating waves are not reflected, a perfect replica 
will be obtained within the wave optics description too 
("perfect lensing") [5]. 



1 



As we have mentioned, in the non-ideal NIM with 
e/x = 1 case considered by Haldane [9] the resonance 
is at a finite wave vector and the repUca is distorted. 
What specifically occurs in the replica image in this case 
has not been discussed. A way to avoid such a discus- 
sion is to redefine the term "resolution" by considering 
the information that, in principle, can be transferred by 
the waves that have passed through the slab [10]. In this 
sense, the resolution limit can be defined as 2'K/kc, where 
kc is the wave vector of the resonant polariton mode. Be- 
low we shall follow a similar definition. 

Consider a slab of a homogeneous negative index ma- 
terial with cjjL = 1, which is parallel to xy plane and 
occupies the region Q < z < I, with one or two sur- 
face layers with thickness (es) s\ and S2- The layers 
also have negative permittivity and permeability with 
eijjLi = e2/i2 = 1- We shall consider below an incident 
p-polarized plane wave such that the magnetic field is 
H = i?o(0, e'('=-'^^+'=-^-'-'*), 0), where and k^ are related 
by Maxwell equations as k^+k^ = ejiu"^ /c^ = w^/c^. Sim- 
ilar results are obtained for s-polarization by replacing e 
by /i in final formulas. The boundary conditions in the 
considered case suggest that both Hy and (dzHy) /e are 
continuous at the interfaces, where e is the corresponding 
dielectric constant. On the "image side" of the slab the 
only nonzero component of magnetic field can be written 
as Hy = ii7oe*'=-(^-')e^('=-^-'^*) (Z ^ / -|- S2 in the case of 
a surface layer), where t is the transmission coefficient, 
which is of main interest to us. In particular, we are 
interested in the transmission coefficient for evanescent 
waves with k^ = in, k = ^k^. — up- jt? > 0, with their 
field Hy decaying away from the slab. 

For a NIM slab with e = — 1 -|- 5, we have 



to\l,d) = e-^'- 



2{l-6) 



sinh kI. 
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If the NIM were ideal {6 = 0), then to{l,0) = expnl, 
and the above mentioned amplification of all evanescent 
waves takes place, which exactly compensates for their 
amplitude decay in the vacuum [5]. In slightly non-ideal 
NIM {S ^ 0) this is already not so, the amplification 

exists only for the waves with < k < Kc 
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where A = 27rc/a; is the radiation wavelength. This es- 
timate is accurate when ^ 1. The pole in the trans- 
mission coefficient to corresponds to a surface polariton 
mode (zero in ^). In the ideal NIM slab this pole corre- 
sponds to Kc = oo, while for slightly non- ideal NIM it is 
at the finite Kc given above [9] . This result tells us imme- 
diately that for any practical slightly non-ideal NIM slab 



the subwavclength imaging is a near- field effect, where 
one basically needs to use very thin slabs comparable 
to the wavelength, Z A, to see any improvement over 
standard focusing [10]. 

Let us now proceed with the study of the effect of the 
surface layers: 

(A) For one surface layer (— si < z < 0) with ei = 
— 1 -I- ?7 and ei/xi = 1, we have the exact answer: 
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sinh kI 



(3) 



(Al) In the case of ideal bulk NIM slab (e = —1), the 
transmission coefficient factorizes: tis = e'^^to{si,r]) [see 
Eq. (1)]. The first factor is associated with the ideal 
NIM slab, while the second one corresponds to the non- 
ideal surface layer. A thin surface layer si <^ X limits 
the resolution as Ax « 27rsi/ In which is much better 
that in the case of the non- ideal bulk (as a whole), cf. 
Eq. (2). 

(A2) In the case of slightly non-ideal bulk NIM the 
transmission coefficient is such that 



tis =to\h6)to\suv) 



5r^{2-6){2-n) 
4(l-^)(l-r?) 



sinh KSi sinh kI. 
(4) 



Now the contributions of bulk and surface non-idealities 
do not factorize, and we find a qualitatively different sit- 
uation. The third term signifies the "hybridization" be- 
tween two polariton solutions that one would have for free 
bulk and surface layers in vacuum. Interestingly, when 
r]5 <Q the polariton pole in the transmission coefficient 
vanishes, because there is a gap in the polariton excita- 
tions at a particular radiation frequency (where €ii = l). 
When the signs of the parameters coincide, r/J > 0, we do 
have a polariton pole. Yet, in both cases the evanescent 
signal is amplified, and the resolution is about the same, 
irrespective of whether there is the polariton resonance 
or not, since the response of the systems remains soft. 

(B) NIM slab with two surface layers. In this general 
case the full solution can also be obtained, but it is a 
bit lengthy and not very instructive. We shall look at 
the most interesting case of ideal NIM slab (e = — 1) and 
surface layers with e\ = —1 + r] and e2 = — 1 -|- ^ where 



t2s = e-^' 
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Again, we may have both regimes, with and without the 
polariton resonance. In particular, if r]9 < 0, the polari- 
ton pole goes off the real axis, but the "amplification" 
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FIG. 1. (a) Inverse transmission coefficient for evanes- 
cent waves (wavenumber k > 0) for a system with e — —1 + 5 
for the bulk slab, ei = —1 + rj, and e2 = —1 + 9 for surface 
layers with si/l = 0.02 and S2/I = 0.01. Curve 1: = 0.05, 
9 — 0.01 with polariton resonance, curve 2: 77 = 0.05, 
9 = —0.01 without the polariton resonance (t"^ / 0), in both 
cases the NIM slab is ideal, 5 = 0. Curve 3: rj — 0.05, 
9 = 0.01, and S = 0.02 (all parts are non-ideal), (b) The 
evanescent field Hy in the replica plane {z — 21) versus the 
field in the source = 0). The evanescent waves with « be- 
low critical value are almost perfectly restored, wave field in 
the region near the polariton resonance is strongly distorted, 
and harmonics with larger k are lost from the replica. 

still exists, Fig. lb. Although the layers are taken as 
very thin, Si/l ^ 0.01, i — 1,2, the suppression becomes 
large due to the cross-terms depending on the thickness 
of the NIM slab I. It is sufScient to consider very small 
deviations from ideal NIM jryj, \6\ ^ 1, when 

t2s ~ e"''('+"i+"'^) (1 - 6'7/e'"(2'+^i+'^2) sinhKSi sinh^sz) , 

(6) 

and we see that the resolution limit is given by 

Ax^ , ^ (7) 

Here, the logarithm is the sum of In / , which is simi- 

lar to the logarithm in Eq.(2) for slightly non-ideal slab, 
and another logarithm of a large argument. In -^j==- The 
two logarithms are expected to be comparable, so in prac- 
tice the situation is not much better than in the case of 
just a slightly non-ideal NIM slab without any surface 
inhomogeneity, cf. Eq. (2). 

The same results as for the surface inhomogeneity 
above can also be obtained within a more general frame- 
work of accounting for a spatial dispersion and the so- 
called "additional boundary conditions" [3]. The effect 



of a spatial dispersion has been mentioned before as one 
of the limiting factors for the perfect lensing, in the sense 
that the subwavelength details in the image will be lim- 
ited by the characteristic length of the microscopic struc- 
ture of the metamaterial acting as a lens, i.e. the spacing 
between the active elements, like metallic split-ring res- 
onators [10]. Our result means that the effects of spatial 
dispersion are much more important than it was recog- 
nized before. 

Prospects for new optical devices, couplers, modula- 
tors, etc. have been driving a strong interest in negative 
index metamaterials. The metamaterials with metallic 
split ring resonators have been used in the first demon- 
stration of negative refraction [4] , and recently their mag- 
netic response frequencies have been extended to 100 
THz (A = 3 /xm) [11] driven by desire to make them 
work in near- visible and visible frequency range. There 
is an ongoing effort in the area of sub-wavelength imag- 
ing with metallic films [12]. Metallic metamaterials ob- 
viously will struggle with losses at optical frequencies. 
This makes nonmetallic systems such as photonic crys- 
tals rather attractive. Photonic crystals made of dielec- 
tric materials have been known to behave as NIM in a 
certain frequency range [13-16], and this was confirmed 
experimentally at microwave [17] and near- infrared fre- 
quencies (A = 1.5—1.6 /im) [18]. It is known that applica- 
tions in this area will face major limiting factors like the 
deviations of bulk parameters from the perfect lens con- 
dition, losses, and a spatial dispersion. In this paper we 
have shown that minute surface inhomogeneities can also 
degrade the subwavelength imaging rather drastically. 
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